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, Matter inhomogeneities along the line of sight deflect the cosmic microwave background (CMB) 

Q\ • photons originating at the last scattering surface at redshift z ~ 1100. These distortions modify the 

0^ ' pattern of CMB polarization. We identify specific combinations of Stokes Q and U parameters that 

correspond to spin 0,±2 variables and can be used to reconstruct the projected matter density. We 
' compute the expected signal to noise as a function of detector sensitivity and angular resolution. 

(£) , With Planck satellite the detection would be at a few a level. Several times better detector sensitivity 

■ would be needed to measure the projected dark matter power spectrum over a wider range of 

scales, which could provide an independent confirmation of the projected matter power spectrum 
as measured from other methods. 



ON 

6 
> 
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Cosmic microwave background (CMB) is believed to originate from the epoch when protons and electrons recom- 
bined into neutral hydrogen making the universe transparent, which happened around redshift z ~ 1100. Photons 
travelling along the line of sight towards the observer were deflected by intervening dark matter distribution via the 
' gravitational lensing effect. This process conserves two quantities: surface brightness and polarization. Conservation 
of surface brightness implies that in the absence of any fluctuations none can be generated by gravitational lensing. 
However, if there are fluctuations present then gravitational lensing can either smooth these out on large scales 
or generate new fluctuations on small scales [||. The effect of these distortions on the CMB anisotropies has been 
thoroughly investigated, both on the CMB power spectrum HD as well as on the induced nongaussian signatures [ffl . 
These studies have shown that the lensing effect is small, but detectable with future CMB experiments, and would 
provide important information on the distribution of dark matter up to z ~ 1100. This method possibly measures 
clustering amplitude at higher redshifts and larger scales than any other method and would be specially valuable in 
breaking the degeneracies between cosmological parameters present when one only uses the CMB data. 

Given the potential importance of measuring the dark matter power spectrum on large scales and/or high redshifts, 
it is worth investigating other methods than can provide similar information. This is particularly important because 
any such measurement will only be statistical in nature and could be susceptible to various systematic effects. One 
such method is gravitational lensing on the CMB polarization. Polarization is also conserved by gravitational lensing 
so that the latter only moves the polarization tensor from one point to another conserving its amplitude and direction 
in the mapping. Since it is widely expected that CMB polarization will be detected with the future CMB experiments 
it is worth asking what kind of information can be learned from investigating the lensing effect on polarization and 
how can it be extracted. 

The effect of lensing on the power spectrum of CMB polarization has already been explored in ||. It has been 
shown there that lensing smoothes the acoustic peaks of E polarization and also generates a B-type polarization. 
While this induced B polarization is small, it provides a fundamental limit to the level of B polarization from other 
sources that can be detected on large scales. For example, it limits the detectable tensor to scalar ratios to be above 
10~ 4 , although this is not a significant limitation for the near future. 

In this paper we explore in more detail the prospects of directly detecting gravitational lensing effect on the CMB 
polarization through its nongaussian signatures. We show that there is a combination of derivatives of polarization 
that provide a local estimator of shear and convergence. These, while noisy, can be used to provide an estimator of 
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projected dark matter power spectrum, which by averaging over all the sky may yield a detectable signal. We compute 
the expected signal to noise as a function of detector sensitivity and apply it to the future satellite experiments. 

II. CONVERGENCE FROM THE CMB POLARIZATION 

In this section we develop a method to extract the gravitational lensing effect on polarization in the ideal case 
without detector noise and beam smoothing. Next section will deal with these complications. A radiation field 
can be fully described in terms of four Stokes parameters || or the specific intensity tensor Iy. For the CMB the 
Stokes parameters can be chosen as the temperature anisotropy T proportional to the sum of intensities along two 
perpendicular directions in the skyQ, the Stokes parameter Q, defined as the difference between the intensities along 
the two axis and the Stokes parameter U, defined by the difference between intensities along the two diagonals. The 
fourth Stokes parameter V describes circular polarization which is not generated via the Thomson scattering believed 
to be the only generating mechanism for polarization, so we will ignore it in the following. 

The Stokes parameters Q and U are two fields in the sky, which we assume to be Gaussian random fields that are 
generated from a common scalar potential (we ignore any initial -B-type plarization, which, even if present, would 
likely be important only on large scales). Because gravitational lensing conserves polarization their observed values 
in the 8 direction Q(8) and U(0) are related to their values at the recombination Q, U, deflected by an angle 88 

Q{8) = Q{8 + 58) 

U{8) = U{8 + S8). (1) 

To extract the lensing information contained in these fields we consider their spatial derivatives. These can be written 
in the weak lensing regime as 

Qa{8) = {5 ab + <5> ab )Q b (8 + 58) 

U a (8) = (5 ab + $ ab )U b (8 + 88), (2) 

where <S> ab — with a,b = x,y are the shear tensor components describing the lensing effect. 

Just like in the temperature case we form quadratic combinations of derivatives of the Stokes parameters and 
express them in terms of unlensed variables and the components of the shear tensor. We wish to begin with a quantity 
that does not depend on the coordinate frame, which allows us to simply relate the unlensed and the lensed quantities. 
Such a quadratic quantity is readily available and is given by Q 2 + U 2 . Then, in the lowest order, we have 

S P = [Ql + Ql + U 2 + U 2 ] (8) 

= (1 + *xx + ®yy)Sp + ($xx - ®yy)Qp + ^xyUp 

Q P ee [Ql -Q 2 y + U 2 - U 2 ] (8) 

= (1 + + ®yy)Qp + - %y)Sp 

Up = [2Q x Q y + 2U x U y ] {8) 

= (1 + ®xx + ®yy)Up + 2<f> xy S P , (3) 

where subscripts x and y stand for respective derivatives in the real space. Equation (||) shows that the measured 
Sp, Qp and Up are products of the shear tensor and derivatives of the unlensed CMB polarization field. Thus the 
power spectrum of Sp, Qp and Up will be a convolution of the power in the CMB and that of the projected mass 
density. Convergence and unlensed fields are taken as independent quantities, which is a good approximation since 
most of the lensing power arises from low redshifts, which are not correlated with the last scattering surface. The 
general expression for this convolution is quite involved. Here we will discuss it in the limit of large scales relative to 
the CMB correlation length £ ~ 0.1° where it is not necessary to take the full convolution into account. In the large 
scale limit patches of the sky larger then the correlation length squared are almost independent, which simplifies the 



1 We work in the small scale limit which simplifies the expressions. Since the final results indicate that lensing on polarization 
is not detectable on very large scales this does not impose a significant limitation on the present work. 
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calculations. This large scale limit is sufficient to analyze future data from MAP and Planck satellites and will, in 
any case, provide a lower limit to the attainable signal. 

The reason we have chosen Sp, Qp and Up as variables in consideration is that they have transformation properties 
similar to the Stokes parameters 0- We wish to show that Sp, Qp and Up transform as Stokes parameters during the 
right-handed rotation of the spatial base (x,y,n) by an angle ip around h directed to the observer. Transformation 
of Stokes parameters Q and U under rotation are given in Q and their linear combinations Q HU have values of 
spin equal to =f2, 

Q' ± iU' = e^> (Q ± iU) . (4) 
We can rewrite definitions of polarization variables (equation ^|) in the following manner 

S P = (Q + iU) x (Q - iU) x + (Q + iU) y (Q - iU) v 

Q P = (Q + iU) x {Q - iU) x - (Q + iU) y {Q - iU) y (5) 
U P = {Q + iU) x {Q - iU) y + (Q + iU) y (Q - iU) x . 

When we use equations (Q) and (||) and change variables in derivatives we find the final expression for transformation 
of our quantities, 

S' P \ / 1 

Q' P = cos 2ip sin 2-0 
Up J \0 -sin 2-0 cos 20 

This demonstrates that Sp transforms as a scalar, while Qp and Up transform like respective Stokes parameters Q 
and U. We will show that these correspond to the familiar convergence and shear quantities of gravitational lensing. 

In the large scale limit (compared to the correlation length of the CMB polarization field) and in the absence of 
lensing we can write the average over the ensemble of the CMB realizations 

{Sp)cMB = 0"5P 
{Qp)cMB = 

{U P )cMB - 0, (6) 

where asp is defined as 

0\SP = (Qt)cMB + (Q 2 y )cMB + (U x )cMB + {Uy)cMB- (7) 

In terms of the polarization power spectrum Cf p it is given by 

aSP = J S (8) 
In the presence of lensing the average of equation (|3|) in the large scale limit becomes 

(Sp)cmb = (1 — "2k)(Jsp 
(Qp)cMB = —1l\OSP 

(Up)cMB = -2l2<?SP, (9) 

where the shear components 71 and 72 and convergence n are defined as 

K = -(^ XX + ^yy)/2 

71 = ~ $ OT )/2 

72 = -$x„. (10) 

Physical interpretation of equation (0) is the following. Convergence k stretches the images and makes the derivatives 
smaller, so (S)cmb is diminished by a factor proportional to k. Similarly, shear produces anisotropy in the derivatives 
in the same way as in the images of distinct galaxies. To normalize the above expressions we introduce 
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such that 

(Sp)cMB = 2k 
(Q'p)cMB = 27i 

{U' P )cMB = 2 7 2- (12) 

In the following we will only use these variables so we drop the primes from now on. 

Instead of working with rotationally non- invariant quantities Qp and Up we combine them to form scalar field Sp 
and pseudo-scalar Bp. In the Fourier space these have the form 

£p{l) = Qp(l) cos(20 z )+Wp(i)sm(20 z ) 

Bp {I) = Q P (l) sin(2^) - U P {1) cos(20 z ), (13) 

where 0j is the azimuthal angle of the mode I. In the Fourier space convergence and shear are related to each other 
through the relations 

71 (Z) = k(1) cos(2&) 72 (0 = k(«) 8111(20,), (14) 
From equations (p]), (|l|) and (|l4|) it follows that 

{Sp)cmb = 2k 
(Sp)cmb = 2k 

<Bp)cmb - 0. (15) 

Thus the convergence k can be reconstructed in two ways, either from Sp or £p. Moreover, vanishing of Bp on average 
can be helpful in indentifing the cosmological part of the signal or removing some foregrounds or other systematics. 
To reconstruct the convergence power spectrum Cf K we form its estimators 

Cjvw' = 1 [ W (iy W '(i') + W(O*W(0] Siv , (16) 

where W or W" stand for Sp, £p, St or St, with the first two obtained from polarization and the last two from 
temperature Q . One can write the mean value of the estimator as 

{C? W )cmb = 4Cf K + N™ w '. (17) 

Here is a power spectrum of the noise arising from the intrinsic fluctuations in CMB temperature or polarization 

(or both for cross-correlation) . This noise arises from the random nature of the CMB and has to be carefully examined 
to be properly subtracted from the estimated power spectrum Cj vw to get an unbiased convergence spectrum C ; KK . 

The noise power spectra can be computed analytically. First, we consider correlation functions between all the 
quantities Sp, Qp, Up, St, Qt and Ut- They are all defined as combinations of derivatives so using equations (||) 
and ( |ll] ) results in the following expressions for the noise 

N SpSp (6) = (Sp(0)S P (6)) 
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i- [iC% Q f + O 2 + {C u Jf + (C u yy u f + 4(Cg u ) 
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[(C£?) 9 ~ (CT y Q ) 2 ] ■ (18) 



Those for temperature are given in j^] and will not be repeated here. The correlation functions C^ y Y are given in the 



Appendix. Equations ( |l8|) together with equations (A2) contain all the necessary information for noise estimate. In 
figure |l| we show correlation functions for polarization and cross-correlation between temperature and polarization. We 
assume throughout that the underlying model is the concordance model with ilf, = 0.04, Vt m = 0.33, SIa = 0.63, 
h = 0.65 and n s = 1. As discussed in |J other viable models give comparable results for signal to noise. For this 
model the correlation length for polarization is £pp ~ 0.07°, while for cross-correlation £pp ~ 0.1° and for temperature 
£pp ~ 0.15°. Each patch of size £ 2 is roughly independent and decreasing £ by a factor of 2 therefore increases by 4 
number of independent patches. Since on large scales CMB noise behaves roughly as a white noise this reduces the 
noise level by corresponding factor. This implies that if we could observe polarization with perfect resolution it would 
give significantly larger signal to noise than temperature alone. However, since polarization has a weaker signal than 
temperature this advantage becomes possible only for very sensitive detectors and small beams, as discussed below in 
more detail. 

The noise power spectra can be obtained from the real space correlations using the following expressions |l3| 

N S A S B = 2jT J edQ n SaSb(q) j q (70) 

Nf ASB =tt J 6d6 {[N QaQb (6) + N UaUb {6)} J (W) + [N QaQb {9) - N UaUb (6)] J 4 (19)} 
n b a b b=1t J ede {[ N QaQb (5/) + N u A u B (0)] ^[n QaQb {9) - N UaUb {9)} J 4 {19)} 

N s A £ B =2tt J 9d9 N SaQb (9) J 2 {19), (19) 
where A and B stand for T or P. From this one obtains the CMB noise power spectra 

N ? pSp = ^f m {2(Qp) 2 + Z{C pp f + 2{C pp f + {C pp f} J (19) 

n[ f£p = ^r-J m {(c<P) 2 + (c! p ) 2 } Hie) + {{c pp f + 2c pp c pp } Mio) 

Nf pBp =^J m {(°o P ) 2 + (C PP ) 2 } MIS) {(C PP f + 2C PP C PP ) J 4 (19) 

N fr£p = _ JL_ f 0dB hc pp C PP + C PP C PP + C PP C PP \ J 2 (W) 
a sp J ' 

n s tSp = f ede | (c tp )2 + 2{c f P) 2 + 4{ c pp ) 2 } mw) 

&ST&SP J > 
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FIG. 1. The upper panel shows auto-correlation functions for polarization N SpSp (9), N QpQp (6), N UpUp (8) and N SpQp (6). 
The lower panel shows the cross- correlations N StSp (6), N QtQp (6), N UtUp (6) and N StQp (6). Both plots are for the cosmo- 
logical concordance model (see text) and without detector noise. 
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We do not need to consider spectra containing Bt or Bp as these pseudo-scalar quantities do not correlate with 
scalars. Fig. || shows noise power spectra for polarization and cross-correlation. The most visible feature is the white 
noise behavior on large scales, which confirms that derivatives of temperature and Stokes parameters are almost 
uncorrelated on scales above a fraction of a degree. 

Let us consider now the power spectra in the limit of low I. Utilizing of the Bessel functions properties in this limit 
and their orthonormality relations we obtain from equation (|20|) the following limits 



lim N? pSp 
i->o 



lim N £p£p 
Yim.N? p£p 



\imN? TSp 



lim N? t£f 
i^o 1 

lim Nf T£F 



= 47T 



(jpdicrr 

lim Nj 3pBp = - lim N? pSf 


47T- 



(21) 



fiddlier? 



lim7V, BTi3p 
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(22) 



The power spectra presented in the fig. g confirm this large scale behavior. It is worth noting that the noise for S 
variables is two times higher than for £ and B, although normalized variances of Sp, Qp and Up are all the same. 
This is because of the spin nature of Qp and Up. The implication of this is that the reconstructed convergence power 
spectrum will have the noise amplitude two times lower when we use £p instead of Sp. 

When we have quantified the intrinsic CMB noise we can use equation ( |l7|) to obtain the best estimate for C^ K from 
observations. In addition, we also need to determine the error on the estimated power spectrum (7 ; ww . To proceed 
analytically we make an additional assumption that Fourier transforms of quantities Sp, £p, St and &r are Gaussian 
distributed. From simulations presented in || this seems to be a good approximation on large scales, where a single 
long wavelength mode receives contribution from many almost independent structures in real space and the central 
limit theorem makes them almost Gaussian. In general, from St, £t, Sp and £p we can construct 10 estimators of 
Cf K , 3 for the temperature field, 3 for polarization and 4 for cross-correlation. Under these assumptions we obtain 
the elements of the covariance matrices of the polarization and cross-correlation estimators 
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FIG. 2. In the upper panel there are noise power spectra for polarization Nf pSp , Nf p£p , Nf pBp and Nf p£p , yet in 
the lower panel for the cross-correlation Nf TSp , Nf T£p , N BtBp and Nf T£p . They are coupled with correlation functions 
presented in figJl via equations 
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where A and £> stand for T or P again. 

To assess the variance of the overall estimator C* K that is a combination of three estimators in the case of polarization 
and four in the case of cross-correlation (eq. [l7]) we treat our variables St, Sp, St and £p as Gaussian. Using the 
Fisher information matrix we can derive the desired variances 111 ill 
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The inverse of the Fisher matrix gives the covariance matrix. The Fisher matrix is one dimensional and relevant 
derivatives are taken with respect to Cf K . We find the following expression for variance of the Cf K estimator when 
we take into account only polarization information 
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Similarly, using the cross-correlation variables only the variance is 
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where we have dropped the cross-correlation terms C WtWp which are much smaller than the diagonal terms in the 
large scale limit. For an experiment like MAP or Planck we will have N~ SpSp ^> N s tSt ^> jyS T s P anc j same 
for other combinations of St, £t, Sp, £p, as can be seen from the noise spectra in figure |]. Only for a significantly 
more sensitive detector will the three noise spectra become comparable in which case, as argued above, polarization 
noise power spectrum will become smaller than that of the temperature. For this reason we provide separately the 
information from polarization and from temperature polarization cross-correlation. In addition, the cosmic variance 
term proportional to ACf K is also much smaller than the CMB intrinsic noise contribution and can be dropped. 
The main source of errors is then the intrinsic CMB noise. In the large scale limit N Sp£p is more than an order 
of magnitude less than N SpSp , N £p£p or j\f B P B p that implies we can neglect that term in the covariances. These 
approximations lead to a diagonal form of covariance matrices that can be written as 
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The minimum variance combination of these estimators in the large scale limit is for polarization 
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In the limit I — > variances of these overall estimators can be written as 
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Var(4Cr) =— - r(Nf pSp ) 2 (41) 
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These expressions show that if N? tSt < Nf pSp then it will be the cross-correlation that will add most of polarization 
information to C Z KK . This information increase will be small compared to the information obtained from the temper- 
ature alone, but may still provide important independent confirmation. If ]\[StSt ^ jyS P s P p i ar i za ^i on provides as 
much information to Cf K as temperature. In further considerations we will use the complete expressions for respective 
variances without the above approximations. 

Another consequence of the intrinsic CMB fluctuations is that we cannot recover the convergence from the individual 
structures, because the coherence length of CMB is too large. We will show that it is still possible to measure the 
convergence in a statistical sense, by averaging over many Fourier modes to extract the power spectrum. We will 
determine the signal to noise for various detector sensitivities in the following section. 



III. CONVERGENCE FROM OBSERVATIONS 



Before computing the expected signal to noise for future missions we need to include the effect of detector noise 
and angular resolution on extracting the convergence power spectrum. Both of these have been discussed in f| and 
the expressions derived there can be applied to polarization as well. 

To quantify the influence of the detector beam on measured Stokes parameters one introduces a filter function F{9), 
which describes the detector beam filtering of the data. The Stokes parameters are convolved with this function 

X a (0) = J F{0- 0') X a {6')d 2 G' (43) 

where X stands for Stokes parameters Q or U. Using equation (||) and Fourier decomposition equation ( ^3| ) becomes 

X a (9) = J d 2 lF(l)X a (l)e d e + J d 2 l J d 2 qF(\l + q\)X b m ab (qy( l +1> e . (44) 

In the presence of filter F(6) it follows from equations (§) and ( [flf ) 

(Sp(9)) cmb = 7A2 / d 2 lF 2 {l)l 2 C[ p l-^Jd 2 q 2K(q)W P (q) e<* 6 . (45) 

We introduced the window function Wp(q) to describe the effect of the detector beam on the convergence. It has the 
form 

, . U 2 lF{l)F{\l + q\)l 2 Cf' P , , 

W P (q) = J - K — K - HU .- 1 . 46 

/ d 2 l F 2 (l)l 2 Cf p 

We see that Wp (q) — ► 1 for q — * independently of the explicit form of the filter. Thus for large scales we reconstruct 
the convergence power spectrum without beam degradation. More generally, the relation between the observables 
and underlying convergence is 

(Sp(1))cmb - (£p(1))cmb - 2k(1)W p (1) (47) 



Window for Planck in the case of polarization and temperature is shown in fig. ga. Just like in the case of 
temperature Q we remove small scale modes with / > l cut where detector noise exceeds primary signal (see below). 
At / ~ 300 the beam already reduces the sensitivity by 50%. Ten times more sensitive detectors would significantly 
extend the range of sensitivty and would be able to measure polarization spectrum up to approximately / w 2000. 
The window for such a detector is also shown in fig. ||a. 

In the presence of the beam smoothing the equation ( |l7|) is rewritten for polarization in the form 

(CF W ') CM B = AW 2 {l)Cr + N™ w , (48) 
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FIG. 3. In the upper panel square of the window function for a few detector configurations using the optimal l cut is shown. 
In the lower panel Nf p£p is shown as a function of l cu t using the white noise approximation. Detector (A) is for current Planck 
proposal. Increasing the detector sensitivity by another order of magnitude, as could be achieved by post-Planck generation of 
CMB satellites, would bring the noise level close to the ideal case. 
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while for cross-correlation its form is 



(Cr W )cMB = iW T (l)W P (l)Cr + N™ w . (49) 

Note that to asses the noise power spectrum we need to measure the CMB polarization, temperature and cross- 
correlation spectrum first. 

Detector noise affects the results through the increase of the CMB correlation length. The detector noise spectrum 
Wp 1 is taken as a white noise in the considered range and is related to the pixel noise and its solid angle size 
through Wp 1 = Cpixei^pixci and similarly for temperature. Planck for example is expected to have the resolution of 
cfwhm = 0.12° and the noise Wp 1 = (0.02) 2 {pK) 2 . Individual structures on small scales are dominated by noise 
and cannot be observed, hence cannot be used to reconstruct convergence so it is best to remove them. We remove 
them by filtering the data on scale corresponding to the scale where noise and signal power spectra are equal, which 
in Fourier space defines some value for l cut . In fact, for the white noise approximation the spectra have the form 
presented in equation ( pl| ) which for polarization gives 



l 5 dl(Cf p +w P 1 e l ' <f 



N, pp = 2tt^ — ^J- -. (50) 

1 iJPdlCf p ) 2 

This can be used to determine l cut which minimizes the noise, at least on large scales where the approximation is valid. 
In fig. |b are shown a few examples of the detector configurations as a function of l cu t . For non-zero detector noise the 
curves have minima in l cut above which noise increases rapidly as expected. Small l cu t filters out real CMB structure 
and increase the correlation length, thus increasing the level of noise. Once l cu t becomes too large the filtering scale 
is too small and small scale power is dominated by noise. This does not have any information on the lensing signal 
and increase the overall noise level again. The optimal value for l cut is found in between these two regimes and agrees 
well with the value defined where the noise and signal power spectra in CMB are equal. 

Figure |J gives the expected temperature, polarization and cross-correlation noise spectra for Planck as well as the 
expected convergence power spectrum. Noise is always several times larger than the expected signal so extraction of 
Cf K is possible only in the statistical sense, by averaging over many mulipoles. In case of polarization the noise is 
about two orders of magnitude larger than the signal. This means that for Planck one would need to average over 
2 x 10 4 modes to obtain S/N ~ 1. This is further reduced by the effect of the window. A more accurate estimate 
of expected signal to noise can be obtained by combining the information from all multipoles and form the optimal 
estimator Q. This gives 

|) -[/^JOS^^f (51) 



ab \ i Var 



(cr) 



where f s k y is the fraction of the sky covered, here taken to be 0.7. Results are shown in the fig. |[ For Planck the 
signal to noise ratio for polarization reaches only (S/N)pp ~ 1 when we take into account about 400 multipoles and 
does not increase beyond that. So the detection of the convergence by Planck from the polarization measurements is 
not really possible, unless the signal in convergence is significantly higher than assumed here. The case of the cross- 
correlation is more promising because of lower noise in temperature. Here we find it is possible to get (S/N)tp ~ 5 
from / < 400. More promising numbers are found if we increase the detector sensitivity by a factor of several, which 
could only be achieved from the next generation of CMB satellites dedicated to polarization. Improving the sensitivity 
by a factor of 10 makes S/N from polarization comparable to the one from temperature, while the cross-correlation 
S/N exceeds both and would provide most of the information. 



IV. CONCLUSIONS 



We have analysed the prospects of measuring the projected dark matter power spectrum up to z ~ 1100 from 
the distortions induced in the CMB polarization. Any such measurement would be extremely difficult both from 
the temperature and polarization because of the statistical nature of such a detection. On the other hand, it would 
provide a direct detection of dark matter clustering on very large scales and high redshifts, possibly not attainable by 
any other method. For this reason it is important that any such detection from temperature fluctuations, described in 
, is independently confirmed. Such a confirmation is possible with polarization, with which one would be provided 
in total with 10 power spectra of convergence, all of which have to agree with each other. 
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FIG. 4. Noise power spectra for temperature (upper panel), polarization (middle panel) and cross-correlation (lower panel) 
measurements. All contain additional noise introduced by Planck detector. 



13 




FIG. 5. Signal to noise ratio as a function of maximum multipole taken into account for temperature autocorrelation (T-T) , 
polarization autocorrelation (P-P) and cross-correlation (T-P). The detector configurations are the same as in figj^ 
. For T-T we use w^ 1 = 0.002^ and <j FWH m = 0.12°. 
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Because polarization is an order of magnitude smaller than temperature fluctuations its detection will be rather 
limited with the planned CMB experiments including Planck. The expected signal to noise for convergence power 
spectrum for Planck is around unity using polarization information only, compared to 15-25 from temperature. Cross- 
correlation could be more promising and could provide S/N ~ 5 or a factor of two larger or smaller depending on 
the actual amount of power in convergence. The real promise of this method lies in the contemplated post-Planck 
experiments dedicated to polarization. Such an experiment would need to cover a significant fraction of the sky to 
measure the projected power spectrum on large scales. If the sensitivity of detectors is reduced by another order of 
magnitude over Planck then measurement of dark matter power spectrum with polarization becomes very significant 
and would extend the amount of information from temperature alone. It would also allow for many cross-checks of 
the results to reduce the possible systematics. 
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APPENDIX A: CORRELATION FUNCTIONS FORMULAS 



Here we present correlation functions of derivatives of the CMB temperature and polarization fields. Due to the 
rotational invariance of the variables St, Sp, £t, £p it is enough to consider correlations between two directions 
separated by 9 along the x direction. During derivation of the correlation functions we make use of the integral 
representation of Bessel functions and the CMB spectra defined in ]l2fl . Weighted spectra Ci(0) are defined as 

cp = J l ^cr jm, (ai) 

where A and B stand for T or P. Some calculations yields the following correlation functions 

(A2) 



CMB 

= (2tt)- 2 j d 2 l e lt - ecos ^ I 2 cos 2 4>iCl P cos 2 2^ 
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= \[Ci P (e) 3/2C pp (9) + C PP {9) 1/2C PP (9)\ 
Cg Q (9) ee (Q y (0)Q y (9)) CMB 

= (2tt)- 2 J d 2 l e «-««»*« / 2 sin 2 <^Cf # cos 2 2<& 

= / S j2C ' ## i [Jo(W) + 3 /2-MW) + J 4 (^) + 1/2J 6 (Z0)] 
= \[C$ P {9) + 3/20^(0) + C pp {9) + 1/2C PP (0)} 

CMS 

= (2tt)- 2 y d 2 Z e «-««»*i ; 2 sin <£, cos <f>iC pp cos 2 20, 
= 

C™{6) ee {U X (0)U X (9)) CMB 

= (2tt)- 2 J d 2 l e «-«coB*, /2 C0S 2^ C PP sin 2 2 ^ 

= / S |2C,## i [Jo(W) ~ V 2 ^) - -M^) + l/2J 6 (W)} 
= \[C$ P {9) l/2C PP {9) C pp {9) + l/2C pp {9)] 
C™<fi) ee (U y (0)U v (9))cMB 

= (2tt)- 2 J d 2 l e il - 6cos ^ I 2 sin 2 4 >l C pp sin 2 2</> ; 

= / S /2c ^i [Jo(w) + l i 2J *w - Jiil6) - l l 2J «W)\ 

= \\ci P {9) + l/2C PP {9) C pp {9) l/2C pp {9)] 

C™{0) EE (U X (0)U V (9))CMB 

= (2tt)- 2 J d 2 l e ^ cos * I 2 sinfa cos 4>iCf p sin 2 2<f>i 
= 

C T J(9) ee (T x (0)Q x (9))cmb 

= (2tt)- 2 J d 2 l e il - ecos(t " I 2 cos 2 4>iCf p cos2<Pi 
= J M f c fpl_ [MW) + MW)] 

= \[ct P {9)-2Ct P (9)+ct P (9)] 
C^(9) ee (f y (0)Q y (9)) C MB 

= (2tt)- 2 J d 2 l e «-0c«*« J 2 sin 2 0,Cf # cos2^ 
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= (2tt)- 2 j d 2 l e 4Z - 9cos * Z 2 sin <^ cos cos 20, 

= 

C™{6) ee {f x (0)U x (9)) C MB 

= (2tt)- 2 J dH e «-0c«*« / 2 cos 2 ^Cpsin 2^ 
= 

C™(#) = (T y (0)f/ 1/ W)CMB 

= (2tt)- 2 J dH e il - 6c ° s *> l 2 sm 2 frcf p sin2fr 
= 

C™(6) = (^(0)f/,W)cMB 

= (2tt)- 2 /" d 2 Z e 4Z - 9cos<# " / 2 sin ^ cos ^Cpsin 2^ 

- f l £ p c? p \[j m - Mm 

= \[CV{6)-Cr{6)] 

CMB 

= (2tt)- 2 y d 2 Z e 4Z - 9cos * Z 2 cos 2 0,Cpsin20 J cos 2<fc 
= 

= (Qv(0)U v (6))cmb 

= (2tt)- 2 y d 2 Z e 4Z - 9cos * Z 2 sin 2 ^Cpsin 2^; cos 2^ 
= 

^(0) = (Qx(O)^(0)>CMB 

= (2tt)- 2 y <i 2 Z e «-»«»*i / 2 sin ^ cos ^Cf p sin 20; cos 2<fo 

= -\[c! p {e)-cp(9)]. 



17 



